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ABSTRACT 

The dynamical aspects of the phonoriton state in highly-photoexcited semiconductors 
is studied theoretically. The effect of the exciton-exciton interaction and nonbosonic char- 
acter of high-density excitons are taken into account. Using Green's function method and 
within the Random Phase Approximation it is shown that the phonoriton dispersion and 
damping are very sensitive to the exciton density, characterizing the excitation degree of 
semiconductors. 
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1. INTRODUCTION 

It is expected that under an intense electromagnetic radiation near the exciton res- 
onance, the interaction between excitons and phonons can lead to the formation of a 
new kind of elementary excitation in semiconductors called phonoriton [1-13]. In treat- 
ing the problem all the authors in [1-13] have considered excitons as ideal bosons and in 
some cases, have neglected the exciton-exciton interaction, although as mentioned therein, 
dealing with phonoriton problem, in fact we have matter with high-excitation phenom- 
ena, where these assumptions are no more valid. It is well known that in highly-excited 
semiconductors the so-called giant polariton states are drastically different from that of 
the usual Hopfield polariton. In such a condition when the number of excitons is very 
large, the excitons cannot be considered as ideal bosons, and the interaction between them 
becomes considerable and cannot be neglected. The nonbosonic character of excitons and 
the exciton-exciton interaction can lead to many interesting effects in semiconductors 
[14-16]. As a mixture of excitons, photons and phonons, the phonoriton spectrum and 
damping must depend on the degree of excitation. However, up to now there has been no 
papers concerning this aspect of the phonoriton problem. 

The purpose of our work is to develop an approach to study the dynamical aspects 
of the phonoriton reconstruction in highly-excited semiconductors. For simplicity, we 
shall use the simple two-level model recently developed by Wang and Birman [11] to 
consider both the above-mentioned effects on the phonoriton state. As will be seen the 
phonoriton spectrum and damping are very sensitive to the exciton density characterizing 
the excitation degree of semiconductors. 

We use the unit system with Ti = c = 1. 

2. THE EFFECTIVE HAMILTONIAN 

Let us consider a semiconductor with direct band gap illuminated by high-intensity 
pump field resonant to the lowest IS exciton state. In such condition the occupation 
number for the exciton in mode k with frequency equal to that of the incident field, 
will be very large (mode k is macroscopically occupied). The creation and annihilation 
operators of the exciton a k , ctk satisfy the following exact commutation relations [16] 



[a k ,a k >] = al,a%, =0 
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+ 
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+ ip*(p + ak)ip(p + ak')h^ +p h k+p } , (1) 

where e+ (e p ), hp (h p ) are the creation (annihilation) operators for the electron and hole 
with momentum p , respectively. ip(p) is the Fourier transform of the exciton envelope 
function for electron-hole relative motion, a = m e / (m e + rrih) = 1 — [3, m e and rrih 
are the effective masses of the electron and hole, respectively. V is the volume of the 
crystal. From (1) it is clear that the exciton can be considered as ideal boson only at low 
excitation, when the second and third terms on the right-hand side of (1) are neglected. 

In addition to the Hamiltonian (1) in Ref.[ll] we must also include the term describing 
the nonlinear exciton-exciton interaction. So our initial Hamiltonian has the following 
form: 

H = H ex ~~< + H ex ~ ph , (2) 
H ex ^ = uj e k x ala k + ujlblb k + g k (alb k + bla k ) + f k alala k a k , (3) 
H ex ~ ph = J2n p c;c p + J2zM(p-q)a;a q (c p . q ~cl p ) , (4) 

V P,q 

where a k (a k ), b k (b k ), c k (c k ) are creation (annihilation) operators for nonbosonic excitons, 
photons and phonons, respectively; uj™ is the energy of the exciton with momentum 
k\io e k x — E t + k 2 /2m ex (E t is the exciton energy at k — 0, m ex is the effective mass of 
exciton), uP k is the energy of the photon , Q p is the energy of the optical phonon, g k is 
coefficient of exciton-photon interaction, M(p — q) is the matrix element of exciton-phonon 
interaction, f k denotes the constant of exciton-exciton interaction. Note that unlike [11] 
a k and a k stand here for nonbosonic excitons and their commutators satisfy (1). Due to 
this nonbosonity of the excitons H ex_1 in (2) can not be diagonalized by the Bogoliubov 
transformation like it was done in [11-13]. To overcome this difficulty we introduce the 
following Green's functions 

Gi(M) = 0(t) < [a k (t),at(o)] > = « a k (t),at(o) » , 

G 2 (k,t) = «b k (t),a+(o) » , (5) 

where the symbol < . . . > denotes the thermodynamical average. 

The equations of motion for the Green's functions (5) can be derived in standard way. 
In order to make the obtained system of equation closed we use RPA-type factorization 
[17] which seems to be reasonable for high density systems where the fluctuation of the 
exciton densities can be neglected. As a result we obtain the closed system of equation 
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for Green's functions in energy representation 

uf -u e k x (l-p k )-2f k n ex (l-p k )]G 1 (k)-g k (l-p k )G 2 (k) = ^~Pk), 

9kGi(k) - (to p k ° l - LdJ)G 2 (k) = , (6) 

where n ex is the exciton density in the semiconductors, and 

Pk = J2{MP- Pk)\ 2 n e (k-p) + \ip(-p + ak)\ 2 n h (k-p)} , (7) 
v 

n e (k),nh(k) are the electron and hole densities, respectively. 

Solving (6) we obtain the dispersion of the polariton i. e. dressed exciton as a pole of 
the Green's function in the following form 

1 

2 



~ ^{M x + 2n e J k )(l-p k )-u;l} 2 + 4(l-p k )gl} 12 , (8) 

i = 1,2. 

In order to study the phonoriton state it is convenient for us to work with bosons 
rather than nonbosonic particles. This idea can be realized if we construct the following 
effective Hamiltonian for the exciton-photon part in term of boson-excitons a k so that 
this Hamiltonian gives us the same polariton dispersion (8) 

HTfT = u e k~ a ta k + uV>tb h + g k (ath + b+a k ). (9) 

Indeed, if we chose u k x and g k , as follows 

£f = (ul x + f k 2n ex )(l -p k ), 



9k = g k \fl-pk, (10) 

we can diagonalize the effective Hamiltonian (9) by a canonical Bogoliubov transformation 
to get it in the form 

KU 1 ' = Y,<kBtkB*k, (11) 

k,i 

where ouf° k l is exactly defined by (8). Bf k is the creation operator for the polariton in 
i-branch and with momentum k, and expressed in terms of photon (b k ) and boson-exciton 
(fife), as follows 

B ik = Ui(k)b k + Vi(k)a k , (12) 



where 



U i( k ) = vol* -y V i( k )> 



pol 

u ik - u k 

V i(k) = T, pol 1 7 n 2 . -211/2' 

IKfe -"I) +9k} y 

Substituting (11) and (12) into (2), and following the same arguments as in [11] we 
obtain the final effective Hamiltonian for the anti-Stokes scattering of the incident macro- 
scopically filled polariton mode k Q of lower polariton branch. 

H AS = £(< + * 7 rH +5 p + E(^ + *7?)s + c P + 

V V 

+ £ tM(p - k )(BZB ko c p _ ko - B£B pC ;_ k J, (14) 

where 

, .pol — pol 
Up = ^2,p > 

M(p-k a ) = M(p-k )v* 2 (p)v 2 (k ). (15) 

Note that the Hamiltonian (14) formally looks like Hamiltonian (12) in [11], however, 
there is an essential difference between them. While the dispersion of the polariton uj v v ° 1 
in [11] is fixed in respect to external excitation, our is function of exciton density 
characterizing the degree of excitation. This is a dynamical consequence of the nonlinear 
exciton-exciton interaction and the nonbosonity of high density excitons. Moreover, we 
have also included lifetime effects introducing r )^° l and 7^ - the damping of polariton and 
phonon, respectively. 

3. PHONORITON SPECTRUM AND DAMPING 

In this section we shall use the effective Hamiltonian derived in the previous section 
to study phonoriton states in highly excited semiconductors. 

Starting now from (14) one can follow [11] replacing Bk a ,B ko by their c-number ex- 
pectation values < B ko > and < B^ o > to get phonoriton dispersion. However, as shown 
in our previous work [13] such a semiclassical approximation may omit many interesting 
properties of phonoritons. So in this paper we shall use the Green's function method [17], 
which is more reasonable and consistent for us in the framework of the RPA. Thus, let us 
consider the following Green's functions 

F 1 (p,t) = «B p (t),B;(0)» , 



F 2 (p,t) = «B ko (t)c p - ko (t),Bt(0)». (16) 

The equations of motion for the Green's functions (16) can be derived in standard way. 
It is easy to check that in the framework of the RPA the Fourier images of the Green's 
functions (16) are the solutions of the following closed system of equations 

' -ctf-i^^M-iMip-kjFfau) = ± 

^ l P-k -illt ko )F 2 {j> ) uj) 



iM(p - ferries VFi(p, + u$° l - 



0, (17) 



where we have approximately replaced the occupation number N ko =< B^ o B ko > of the 
polariton at mode k Q by N ex = V.n ex due to the exciton-like nature of the polariton mode 
k D near the exciton resonance. 

From the solutions of the system of equations (17) we extract the poles of the functions 
Fi(p, uj) and F 2 (p, cu), the real and imaginary parts of which, as usual, give us the spectrum 
and damping of the new quasiparticle - the phonoriton. 

For the spectrum of phonoriton we obtain the following final expression 
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and for the damping of phonoriton we have 

r«(p) : 



where 
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(20) 
(21) 
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- {[« + 2f k n ex ){\ - p k ) - ujJ] 2 + AgUl - p k )} ] 



(23) 



p k \ir.a B .n ex [16], where A is some constant parameter. The formulas (18) and (19) are 
the final expressions for the dispersion and damping of the phonoriton. As we see, these 
analytical expressions depend explicitly on the exciton density that characterises the level 
of excitation. From (18) we have a new gap of phonoriton 

A(p) = |^[ J D(p) 2 + 4M 2 (p-A; )yn e:c -7 2 (p)" 



+ 



1 

+ 2 



D{ P y + m\ V - k )Vn ex - 7 2 (p)l 2 + AD{ P y^{p) (24) 



It is clear from (18)-(19) that the phonoriton dispersion and damping depend on 
exciton density not only via AM 2 (p — k )Vn ex as in [10-12], but also through the polariton 
dispersion (23). If we put fk = and p k = we reobtain the results of [11]. 

To see the sensibility of phonoriton states to the exciton density, in Fig.l, 2, and 3 
we have plotted the phonoriton dispersion, its gap and damping in CdS for some exciton 
densities, respectively. The damping of polaritons can be calculated via the damping 
of excitons and photon (see for example [18]), and / = l3irE B a 3 B /3V [19], where a B is 
the effective Bohr radius of the exciton and E B is its binding energy. The data used 
for numerical calculation are adopted from [11]: E t = 2553mev, g Q = 50mev, a B = 
2.8 x l(r 7 cm, E B = 27.8mev, k Q = 3.98 • 10 5 cnT 1 ,m ex = 0.98m o ,Y x = O.lmeV, Y h = 
0.05meV. Note that the higher exciton density the lower two phonoriton branches go 
down. With increasing exciton density not only the phonoriton gap becomes larger as 
in [11], but also the phonoriton resonance region shifts to higher values of momentum. 
This situation may affect experimental detection of phonoriton and must be taken into 
account. 



4. CONCLUSION 

In this paper we study for the first time the dynamical aspects of the phonoriton states 
in highly excited semiconductors. Unlike the previous papers on phonoritons considering 
the polariton as fixed in respect to external excitation before it comes to interaction with 
phonon, we have shown that the polariton itself depends on excitation. The intrinsic non- 
linear interactions to high-density systems such as nonlinear exciton-exciton interaction 
and nonbosonic nature of excitons are taken into account. As consequence, these nonlin- 
ear interactions lead to dynamical reconstruction of phonoriton states. At fixed k -mode 
of photon excitation, the higher exciton density (the higher excitation degree) the deeper 
the phonoriton resonance region goes down, and the larger values of momentum k the 
latter shifts to. Note that if we neglect these interactions, i.e. f k — * and p k — > all 
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our results return to the previous ones of [11]. The sensibility of the phonoriton spectrum 
and damping upon the degree of external excitation usually forgotten up to now should 
be taken into account, especially in experimental detection of phonoritons. And we hope 
that our specific predictions will stimulate further experiments. 
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FIGURE CAPTIONS 



Fig.l The dispersion of the anti-Stokes phonoriton at some exciton densities n ex = 1-10 
cm -3 (thick line), 1.5 • 10 17 cm -3 (thin line), 2 • 10 17 cm" 3 (dotted line). 



Fig. 2 The function A(p) at some exciton densities n ex = 1-10 17 cm 3 (thick line), 1.5-10 
cm" 3 (thin line), 2 • 10 17 cm" 3 (dotted line). 



Fig.l The damping of the anti-Stokes phonoriton at some exciton densities n ex = 1-10 
cm" 3 (thick line), 1.5 • 10 17 cm" 3 (thin line), 2 • 10 17 cm" 3 (dotted line). 
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